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ABSTRACT 

A  generalized  waterjet  propulsion   system  in   a    sub- 
cavitating  hydrofoil   craft   is    considered   for  design   and 
performance   estimation.      Independent   system  parameters   of 
jet  velocity   ratio,    inlet  velocity   ratio,    and  nacelle   inlet 
diameter   to  maximum  diameter   ratio    are    varied   to   search    for 
the    optimum   system.      The   ODtimum  system   is   defined   as   the 
minimum  total  propulsion   system  weight.      The    optimization 
scheme    utilizes   a   directed   search  without  the    calculation 
of   derivatives,    and  was    chosen    for  its    simplicity,    versa- 
tility,   and   rapidity.      The    ducting   system  is   divided   into 
components    for  head   loss   calculations.      For   a  given   design, 
total  head   losses   in   the    duct    are    computed  by   means   of 
experimental   data   and  empirical  equations,    enabling  pump, 
reduction   gear,    and  prime   mover  desicrn  to  be    completed. 
Results    for   a  sample    craft   are    included  which  indicate   that 
hydrofoils   should  be    designed   about  the   gas    turbine,    due    to 
the   discrete  power   levels   that  are    available.      Suboptimizing 
the  nacelle   desion  may   lead   to  an   overall    less   optimum 
system  design.      A  general    lack   of   information  was  noted   on 
nacelle,    cascade    corner,    and  three    dimensional   diffuser 
design   and  performance   prediction.      A  FORTRAN   computer 
listing   and   flow   charts    are    included. 
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CHAPTER  1  -  INTRODUCTION 
The  use  of  waterjets  for  the  propulsion  of  marine 
vehicles  at  conventional  ship  speeds  has  not  been  consid- 
ered comparable  to  that  of  propellers.   Although  the  con- 
cept is  quite  old,  the  disadvantages  of  waterjets  have  not 
been  overcome  and  the  advantages  not  realized.   At  the 
speeds  considered,  the  low  system  efficiency,  large  system 
weight  and  relatively  sophisticated  technology  required 
were  sufficient  to  virtually  exclude  waterjets  from  com- 
petition with  propellers.   With  the  advent  of  high  speed 
craft  however,  this  means  of  propulsion  is  receiving  a 
renewed  interest  as  a  practicable  system.   Among  the  pos- 
sible advantages  that  may  be  realized  are  (refs.  1,2): 

a.  Reduced  noise  and  vibration 

b.  Fewer  system  components 

c.  Elimination  of  underwater  appendages 

d.  Elimination  of  complex  transmission  machinery 

e.  Limited  draft  applications 

f.  Simplification  of  steering  devices 

For  hydrofoil   craft,    the  waterjet   looks  particularly 

attractive.      It  neatly   avoids   the  problem  of  conventional 

propellers    in  power  transmission    (refs.    3,4).      For  this 

reason,   waterjet  propulsion   development   to-date   has  been 

primarily   concerned  with  hydrofoil   applications.      In 

addition,    new   concepts    in   ships,    such   as   surface   effect 
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ships  and  captured  air  bubble  craft,  are  now  including  con- 
sideration of  waterjets  as  a  possible  propulsion  system. 
With  this  surge  of  interest,  the  need  to  go  beyond  the 
practical  "Will  it  work?"  engineering  design  method  has 
become  paramount.   The  question  has  become,  "What  is  the 
best  waterjet  propulsion  system  for  the  craft  under  consid- 
eration?".  Answering  this  question  requires  a  strictly 
formulated  methodology  previously  not  required  in  make-it- 
work,  technological  levels.   No  longer  is  it  satisfactory 
to  find  some  point,  any  point,  above  the  minimum  acceptable 
level.   The  scope  has  been  enlarged  to  pinpointing  the 
optimum  within  the  entire  domain  of  feasibility.   The 
immensity  of  this  task  coupled  with  its  relative  lack  of 
need  has  hampered  the  implementation  of  a  methodology. 
Hydrofoils  have  demonstrated  the  need,  while  computers 
provide  the  invaluable  aid  for  the  required  computations; 
therefore,  an  optimization  methodology  has  been  developed. 
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CHAPTER  2  -  DESIGN  METHODOLOGY 

To  initialize  the  design  method,  a  clear  definition 

of  the  desired  objective  must  be  given.   In  the  present 

instance,  the  definition  is  the  answer  to  the  question, 

"What  is  meant  by  an  optimal  water jet  propulsion  system?". 

The  choice  of  an  optimum  must  adequately  describe  the 

system  without  bias.   It  is  readily  apparent  that  whatever 

optimum  is  chosen,  it  will  have  a  strong  effect  on  the 

methodology  employed.   The  current  most  commonly  used 

descriptor  is  the  propulsive  efficiency  which  is  defined 

as  the  ratio  of  the  power  required  to  the  power  output  of 

the  plant  (ref.  5).   Efficiencies  suffer  from  differences 

in  definition  and  thus  vary  in  their  relative  importance  to 

system  optimization.   Minimum  power  required  to  propel  the 

craft  at  the  desired  speeds  is  sometimes  utilized  as  an 

optimization  objective  with  the  advantage  of  being  readily 

calculable  and  understood.   A  disadvantage  is  that  a  gas 

turbine,  which  is  a  commonly  used  prime  mover  on  account  of 

its  low  weight  to  power  ratio,  operates  most  efficiently  at 

its  maximum  power  level.   Hence,  there  exists  a  natural 

trade-off  between  the  most  efficient  propulsion  and  the 

margin  required  for  growth  in  the  ensuing  years  of  the 

craft's  life.   Minimum  power  calculations  are  certainly 

necessary  but  provide  a  somewhat  questionable  optimum.   In 

an  attempt  to  avoid  vagaries  of  definition  and  usage  and 
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to  remove  doubts  of  sufficiency,  the  goal  of  optimization 
or  the  system  "cost"  may  be  defined  differently.   The 
"cost"  of  any  system  should  reflect  the  disadvantages  of 
that  system  without  weighting  one  disadvantage  over 
another.   Two  disadvantages  of  waterjet  propulsion  systems 
mentioned  previously  are  a  large  system  weight,  in  part 
due  to  the  entrained  water,  and  a  low  overall  system 
efficiency.   Inefficiencies  are,  however,  necessarily 
manifested  in  increased  weight,  be  it  added  required 
fuel,  another  pump  stage,  or  larger  ducting.   Recall,  too, 
that  hydrofoils  are  weight  limited  vehicles.   A  unique 
number  that  reflects  these  realizations  is  a  suitable 
definition  of  the  weights  attributable  to  the  propulsion 
system  (ref.  5).   The  propulsion  system  weight  shall  be 
defined  as  the  sum  of  all  weights  of  the  structural  com- 
ponents, prime  movers,  fuel,  pumps,  gearbox,  ducting,  and 
entrained  water  above  the  mean  waterline.   If,  for  a  given 
craft,  the  displacement  is  considered  constant,  then  minimi- 
zation of  the  propulsion  system  weight  is  equivalent  to 
maximizing  the  craft  payload.   Thus,  a  propulsion  system 
shall  be  defined  as  optimal  if,  in  addition  to  fulfilling  the 
propulsion  requirements,  it  adds  the  minimum  weight  of 
any  propulsion  system  to  the  vehicle. 

The  generalized  waterjet  propulsion  system  for 
hydrofoils  utilized  in  this  study  is  depicted  in  Figure  1. 

It  should  be  noted  that  this  configuration  permits 
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FIGURE  1    GENERAL  DUCTING  CONFIGURATION 
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considerable   variation  by   the    changing   of   a   few   lengths 
and  angles   of   the    ducting.      The  water  enters   the   nacelle 
under  ram  pressure   and   is   ducted  up   the    strut  and  into   the 
pump  where   the  energy  input  is   converted   to  kinetic 
energy   in   the    jet.      Thus   the    thrust   reaction   is   generated, 
providing   the    craft's   propulsion.      In  equation    form, 

T  =    pQ(Vj    cos (6)    -   V0)  (1) 

In  any  optimization  procedure,  the  formulation  of  the 
method  requires  the  knowledge  of  the  parameters  to  be 
varied.   The  number  and  selection  of  these  parameters 
contributes  significantly  to  the  applicability  of  the 
method.   The  parameters  must  be  necessary  and  sufficient 
system  descriptors,  and  it  is  desirable  to  utilize  the 
minimum  number  required.   If  M  operations  are  performed  in 
the  evaluation  of  one  design  point  and  N  independent  para- 
meters describe  the  system,  then  the  total  number  of 
operations  required  increases  as  M  .   The  sufficiency  of  the 
description  of  the  system  is  clearly  a  choice  of  adequacy. 
The  variables  that  most  strongly  affect  the  resulting  system 
should  be  chosen  as  system  parameters  while  other  variables 
that  produce  only  minor  variations  over  a  wide  range  of 
values  on  the  essential  design  should  be  defined,  or  even 
somewhat  arbitrarily  chosen,  to  suit  the  needs  of  the  design. 

Examination  of  equation  (1)  reveals  that  for  a  given 

craft  and  associated  mission,  only  the  flow  rate  and  jet 

velocity  are  unknown.   The  one-dimensional  continuity 
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equation    for   an   incompressible    fluid 

Alv!  =   A2V2   =   Q  (2) 

indicates    that   choice   of   jet   velocity  will   determine    the 

jet   area.      The    first  non-dimensionalized   system  parameter, 

then,    shall  be   V-/V0,    defined  at   the    cruise   point.      Knowing 

the    flow   rate    from   the    jet,    the    same    reasoning  may  be    applied 

to   the    inlet  which   results   in    the    second   system  parameter  of 

V^/VQ,    where  V^    is    defined  in   the    cruise    condition   at   the 

point   of  minimum  area   in   the   nacelle.      The   third   and 

last  system  parameter,    for   reasons   that  will  be   explained 

in   Chapter    3,    shall  be    D^/Dm,    where   Di    is    the   diameter  of 

the   inlet  at  the  point  of  V^ . 

Finally,    the   minimum   amount  of   information  needed  to 

specify   the    craft   and   its  mission  must  be    determined.    . A 

typical   drag   to  craft  speed  curve    is   shown   in   Figure   2. 

The    "hump"    at  the    lower   speeds   represents   the    take-off 

point,    i.e.,   when  the   hull   lifts    free   of  the  water  and  the 

craft  is    foilborne,    and  generally   occurs    at   speeds   that   are 

approximately   50%   of   the   cruise   speed.      The   sizing   of   the 

propulsion   system  then  is   dependent  on  both  the   cruise   and 

take-off   conditions.      The  minimum  information   required  is 

the  following: 

a.  Displacement 

b.  Range 

c.  Prime  mover 

d.  Craft  SDeed,  at  take-off  and  cruise 
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e.  Craft   Drag,    at   take-off  and  cruise 

f.  Depth  of   submergence   of   foil 

g.  Distance   of   strut   from  transom 

h.      Distance   of  pump  exit   from  transom 
i.      Height  of  pump   centerline    above   keel 
Although  the  prime  mover  is   assumed  to  be   a  gas   tur- 
bine,   knowledge   of   the  particular   type   of  gas   turbine   is 
important  in   consideration   of  both    fuel   and  power  re- 
quirements.     It  is   readily   seen   that  the    last   four   input 
requirements   define   the   approximate   configuration  while 
the    first   five    items   describe   the   craft  size   and  mission. 
Knowing   the    required  ingredients    and   the    desired  re- 
sult  leaves   only   the   recipe   to  be   set   forth.      For  the 
traverse   of  an   ideal   fluid  between  two  points   in   a  closed 
conduit,    Bernoulli's   equation  may  be   applied  such   that 

Pi   +    2pV?   +    pgzl  =  ?2   +   IpV2   +    pgz2  (3) 

where    2   indicates   a   station   downstream  of   station   1. 

For  a   real   fluid  with   losses   this  becomes 

?  9 

pl  +  2PVX  +  P9zi  "  P2  +  2^v2  +  pg22  "  Ap12   (4) 
where  Ap,2  =  losses  between  stations  1  and  2. 

Although  this  equation  can  be  applied  only  along  a  stream- 
line and  hence  cannot  be  used  across  the  pump,  it  is 
applicable  to  the  ducting  on  either  side  of  the  pump.   The 
continuity  equation  provides  the  information  that  the 
ducting  velocities  are  determined  by  geometric  configura- 
tion.  Starting  from  a  known  point,  e.g.  the  water  surface, 
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and  following  the  water  through  the  duct  leaves  only  the 
losses  as  unknown.   Estimation  of  the  losses  may  be  accom- 
plished based  on  duct  size  and  shape  (ref .  8) .   This  yields 
the  needed  information  for  the  pump  design. 

Following  the  flow  diagram  of  Figure  3,  and  recalling 
the  previous  discussion,  the  choice  of  the  system  parameters 
by  the  technique  described  in  Appendix  A  determines  the  in- 
let and  outlet  geometry.   Knowledge  of  the  jet  area  determines 
the  flow  rates  at  both  cruise  and  take-off.   Water  is  in- 
gested into  the  nacelle  at  velocity  V^ .   The  parameter  D-;/D 
sets  the  amount  of  diffusion  possible  prior  to  entry  into  the 
strut  elbow.   Determination  of  the  nacelle  length  provides 
both  a  drag  estimate  and  the  head  loss  in  the  nacelle  duct 
as  described  in  Appendix  B.   The  strut  elbow  will  exhibit 
large  losses  due  to  the  low  aspect  ratio  (width/depth)  of 
the  duct  as  it  conforms  to  the  shape  required  for  the  strut 
and  due  to  the  relatively  large  velocities.   The  use  of 
guide  vanes  will  minimize  the  losses  by  increasing  the 
aspect  ratio.   Assuming  the  inlet  and  outlet  areas  are  the 
same,  the  losses  in  the  elbow  may  be  calculated  by  the  method 
described  in  Appendix  C.   The  strut  diffuser  area  ratio  is 
set  to  avoid  cavitation  in  the  downstream  ducting  and  to 
reduce  the  inflow  velocity  to  the  pump.   Sizing  the  internal 
flow  area  permits  an  estimate  of  the  exterior  dimensions 

which  are  sized  up  to  cavitation  limits.   For  more  information, 
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see  Appendix  D.   The  drags  attributable  to  the  strut  and  spray 
are  estimated.   The  newly  calculated  nacelle,  strut,  and 
spray  drags  are  compared  to  the  original  estimates.   If  the 
deviation  is  excessive,  the  drag  estimate  is  revised  and 
the  problem  restarted.   Otherwise,  the  head  losses  in  the 
strut  due  to  friction  and  diffusion  are  calculated.   Entry 
into  the  hull  from  the  strut  is  made  via  a  90°  elbow  for 
which  the  losses  are  calculated  as  before.   It  should  be 
noted  that  this  elbow  will  have  a  more  favorable  aspect 
ratio  as  a  result  of  the  strut  shape  and  will  require  fewer 
guide  vanes.   Coupled  with  the  decreased  velocity,  the 
losses  will  show  a  decrease  from  those  due  to  the  strut 
elbow.   Determination  of  the  ducting  required,  as  described 
in  Appendix  E,  permits  estimation  of  the  ducting  losses 
therein.   Summation  of  all  the  losses  prior  to  the  pump  and 
knowledge  of  the  inflow  velocity  provide  the  cavitation 
information  in  the  form  of  net  positive  suction  head  (NPSII) 
to  the  pump.   Adding  the  nozzle  losses  permits  calculation  of 
the  pump  head  required.   The  method  for  determining  the 
nozzle  losses  may  be  found  in  Appendix  G.   This  is  sufficient 
to  determine  the  size  and  type  of  pump  needed,  the  details 
of  which  may  be  found  in  Appendix  F-   The  specifications  of 
the  prime  mover  are  then  utilized  to  specify  the  gearbox 
parameters  of  gearbox  ratio  and  power  which  in  turn  indi- 
cates the  SFC  required  in  the  cruise  condition.   Since  the 
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amount  of  fuel  required  for  the  endurance  specified  is 
calculable,  the  entire  propulsion  system  is  fully  defined 
in  dimensions,  configuration,  weights,  and  particular  com- 
ponents.  Summing  all  the  component  weights  yields  the 
desired  system  weight.   This  concludes  the  evaluation  of 
a  propulsion  system  weight  for  the  chosen  system  parameters, 
The  optimization  routine  utilizes  the  weight  information, 
selects  a  new  set  of  system  parameter  values  and  the 
cycle  begins  anew.   Ultimately,  the  minimum  weight  system 
is  located  and  the  orocess  completer. 
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CHAPTER  3  -  DISCUSSION  OF  RESULTS 
A  sample  output  is  depicted  in  Figure  6.   The  first 
page  is  predominantly  a  repetition  of  the  input  data  and 
is  printed  out  upon  call  to  the  program.   The  second  and 
third  pages  are  the  resulting  optimized  system  data  that  is 
printed  for  each  configuration  of  pumps  and  gas  turbines 
considered.   Hence  the  designer  can  see  at  a  glance  the 
pertinent  information  required  for  a  preliminary  design  as 
well  as  compare  the  alternative  configurations  that  might 
be  feasible  and  competitive  designs. 

There  are  a  number  of  features  within  the  results  that 
are  not  immediately  obvious.   First  and  foremost  is  the 
problem  of  results  exceeding  a  system  limitation.   For 
example,  having  required  cruise  shaft  horsepower  exceeding 
the  maximum  normal  horsepower  of  the  engine.   Within  the 
program  during  the  optimization,  such  designs  are,  of 
course,  rejected.   For  the  print  out,  it  was  considered  that 
knowledge  of  the  resulting  "optimum"  was  desirable,  how- 
ever erroneous  it  may  be.   One  of  the  forcing  influences  on 
this  decision  was  the  variability  of  the  total  drag  estima- 
tion and  thus  the  flow  rate.   If  the  optimum  system  re- 
sulted in  a  value  that  was  at  a  boundary,  a  slight  change 
in,  say,  flow  rate  might  cause  the  value  to  exceed  the 
boondary  and  hence  result  in  a  rejected  design.   In  other 
words  the  optimum  system  is  not  identically  reproducible. 

The  significance  of  this  fact  will  be  amplified  in  the 
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conclusions  section.   This  lack  of  reproducibility  results 
primarily  from  the  iterations  for  drag  estimation.   Each 
new  estimation  is  compared  to  the  previous  estimation  and 
accepted  if  within  5%  of  the  former  value.   Otherwise  the 
flow  rate  is  adjusted  for  the  new  drag  and  the  nacelle  and 
strut  are  redesigned.   In  the  case  of  the  400  ton  craft 
used,  5%  of  the  total  drag  is  on  the  order  of  4000  pounds 
force  in  the  cruise  condition.   Thus,  starting  from  the 
same  system  parameter  values  but  with  different  initial 
drag  estimates  may  result  in  a  total  drag  differing  by  5% 
in  the  output  with  corresponding  differences  in  the  system. 

Secondly,  the  resulting  design  is  obviously  constrained 
by  the  ducting  configuration  used.   Hence  direct  comparison 
with  existing  craft  is  not  valid  until  the  configuration 
is  modified  to  conform  to  that  of  the  craft. 

At  the  inception  of  this  project,  it  was  believed 

that  the  duct  area  ratio  of  the  strut  was  a  system  parameter 

to  be  varied  independently.   Clearly  a  lower  bound  on  this 

parameter  is  that  amount  of  diffusion  required  to  avoid 

cavitation  in  the  hull  ducting.   It  was  felt,  however, 

that  the  pump,  gear  and  fuel  weight  calculations  might 

reflect  a  lower  overall  weight  if  the  pump  inflow  rate  was 

reduced.   The  results  did  not  justify  this  assumption.   The 

system  unerringly  chose  the  strut  with  the  minimum  required 

diffusion,  which  meant  minimum  weight.   The  principal 

reason  for  this  result  is  that  the  pump  design  is  in  terms 
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of  total  head.   Hence  the  pump  effectively  never  "sees"  the 
amount  of  static  head  presented  to  it.   The  pump  is  affected 
only  by  the  differences  in  ducting  losses  which  are  re- 
flected in  both  the  NPSH  and  pump  head.   Thus  the  net  effect 
is  of  secondary  importance  and  is  not  significant  enough 
to  change  the  area  ratio  of  the  strut.   For  this  reason  the 
area  ratio  of  the  strut  ducting  was  rejected  as  a  system 
parameter  and  is  presently  determined  as  required  by 
cavitation  limitations  or  pump  inlet  area. 

The  initial  nacelle  design  method  chose  the  largest 
inlet  to  maximum  diameter  ratio  permissible  by  exterior 
cavitation.   This  ratio  was  tempered  by  the  amount  of 
diffusion  required  to  avoid  cavitation  in  the  strut  elbow. 
Assuming  the  same  inlet  velocity  ratio,  would  more  diffusion 
than  the  minimum  required  permit  an  overall  lesser  weight 
system?   In  other  words,  does  the  optimum  system  require 
the  optimum  nacelle  or  might  an  off  design  nacelle  be 
used  resulting  in  a  lesser  weight  system?  The  use  of  the 
nacelle  diameter  ratio,  DIDM,  as  a  system  parameter  pro- 
vides the  opportunity  to  investigate  this  question.   The 
current  results  are  somewhat  inconclusive  due  to  the  influence 
of  other  factors.   Tentatively,  it  is  suggested  that  for  some 
designs,  the  use  of  DIDM  as  a  system  parameter  is  justi- 
fied and  results  in  a  better  system  design.   It  is 
recommended  that  investigation  of  this  question  be  pursued 

for  determination  of  what  constitutes  the  basis  of  nacelle 
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design  in  regards  to  system  performance. 

Perhaps  the  most  significant  conclusion  is  that  the 

choice  of  engine  is  extremely  important  in  the  resulting 

optimum  system.   The  first  reason  for  this  is  that  the 

gas  turbine  is  presently  available  only  at  discrete  power 

levels.   Since  fuel  weight  is  large  relative  to  the  other 

system  weights,  the  operation  of  a  gas  turbine  at  a  point 

significantly  off  its  design  point  will  result  in  a  large 

increase  in  fuel  weight  and  thus  system  weight.   The  net 

effect  is  that,  given  the  opportunity  to  avoid  the  penalty 

in  fuel  weight,  the  system  tends  towards  the  maximum  SHP 

permitted  by  the  engine,  and  will  drive  the  jet  velocity 

ratio  as  high  as  possible  to  achieve  this  SHP.   This  also 

results  in  a  low  flow  rate,  low  drag  and  low  propulsive 

coefficient.   If  the  engine  is  slightly  large  for  the  craft, 

then  the  results  will  show  a  high  jet  velocity  ratio,  SHP 

per  engine  at  the  maximum  permissible  level,  which  may  be 

limited  by  take-off  or  cruise  SHP,  and  a  low  propulsive 

coefficient.   If  the  engine  is  so  large  that  there  is  no 

power  limitation,  the  jet  velocity  ratio  will  be  very 

high  and  the  SHP  per  engine  relatively  low  resulting  in  a 

very  low  propulsive  coefficient  due  to  large  total  SHP 

and  a  high  system  weight.   The  effect  of  an  engine  that  is 

barely  adequate  is  a  lowered  jet  velocity  ratio  in  an  effort 

to  reduce  the  SHP  requirements.   The  trade-off  here  is  the 

large  flow  rate  required  to  compensate  for  the  low  jet 
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velocity  ratio  and  the  large  drag  that  results  from 
attempting  to  accomodate  the  flow  rate  in  the  nacelle  and 
strut.   Thus  SHP  is  being  forced  up  and  the  system  may  never 
achieve  an  acceptable  horsepower  level. 

This  power  level  sensitivity  is  so  strong  that  it  is 
abundantly  clear  that  hydrofoils  should  be  designed 
around  a  particular  gas  turbine.   Choice  of  a  desired 
displacement  and  the  use  of  an  adequate  prime  mover  may 
result  in  a  poor  design  and  inefficient  use  of  the  gas 
turbine.   Choice  of  the  gas  turbine  and  permitting  the 
displacement  to  be  chosen  at  the  point  of  best  engine 
operation  will  result  in  a  more  efficient  and  more  effective 
design.   Thus  it  is  recommended  that  another  version  of  this 
program  be  created  that  chooses  the  best  displacement 
given  the  engine.   It  is  also  suggested  that  perhaps  the 
concept  of  hydrofoil  craft  design  be  revised  to  accomodate 
the  fact  that  the  engine  choice  determines  the  best  design. 
Figure  7  shows  the  minimum  weight  ratio  for  one  LM2500 
and  one  pump  with  a  lift  to  drag  ratio  of  11.3.   The 
results  show  the  optimum  weight  ratios  with  and  without 
fuel  do  not  coincide.   The  propulsive  coefficient  is  drawn 
to  show  the  trade-off  between  maximizing  SHP  and  maximizing 
the  propulsive  coefficient. 

The  choice  of  an  engine  also  affects  the  reproducibility 

of  the  optimum  system  design.   If  the  engine  is  smaller 

than  desirable,  VJVO  is  driven  as  high  as  the  SHP  permits. 
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When  the  optimum  system  parameters  have  been  found  and  the 
optimum  system  is  to  be  regenerated,  the  design  process  is 
reinstituted.   During  this  final  design,  iterations  in 
craft  drag  are  again  required  to  determine  the  final  flow 
rate  and  size  of  the  system.   Thus  the  net  result  may  be 
a  system  slightly  different  from  the  one  originally  found 
with  least  weight.   And  in  the  case  of  a  system  already  on 
a  boundary,  the  resulting  system  may  exceed  the  limits 
placed  on  the  optimization  process.   It  is  to  be  expected 
that  the  deviation  will  be  small  but  nonetheless  present. 
The  results  in  such  a  case  should  be  interpreted  as  optimal 
at  the  boundary  and  not  an  unsatisfactory  design.   If  the 
deviation  is  significant  and  the  other  system  values 
reflect  it,  the  interpretation  should  be  that  the  system 
was  unable  to  be  optimized  within  the  limits.   The  most 
common  example  of  this  occurs  when  the  cruise  or  take-off 
SHP  of  the  engine  is  insufficient  to  meet  the  craft  re- 
quirements. 

The  effect  of  range  on  the  system  design  is  shown  in 
Figure  4.   The  LM  2500  is  a  power  limited  engine  for  the 
four  hundred  ton  craft  at  forty-two  knots.   The  FT4C-2, 
however,  is  so  large  that  maximum  power  is  never  attained  in 
this  craft.   The  effect  of  lower  SFC  and  operation  at 
maximum  power  favors  the  LM  2500  at  the  longer  ranges.   As 
it  should  be,  the  weight  ratio  without  fuel  is  essentially 

constant  for  both  engines,  although  the  weight  ratio  for  the 
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LM  2500  is  higher,  presumably  on  account  of  the  power  limita- 
tion.  Similarly,  the  jet  velocity  ratio  remains  constant 
as  does  the  flow  rate,  but  the  inlet  velocity  ratio  tends 
to  increase  with  increasing  range.   The  nacelle  diameter 
ratio,  on  the  other  hand,  shows  no  consistent  pattern. 
The  effect  of  variation  of  displacement  by  10%  is  shown  in 
Figure  5  for  the  FT4C-2  for  one  and  two  pumps. 

The  results  show  a  surprising  sensitivity  to  the 
starting  values  of  the  system  parameters.   When  arbitrarily 
deviated  by  twenty  per  cent  from  the  present  values,  the 
results  were  unpredictable.   The  current  values,  however, 
have  shown  good  stability  with  test  craft  ranging  from 
fifty  to  seven  hundred  fifty  tons  and  from  forty  to  fifty 
knots.   The  reason  for  this  apparent  instability  is  un- 
known.  The  step  sizes  appear  to  be.  acceptable  at  any  size 
above  a  threshold  value.   Some  conditions  will  not  permit 
location  of  the  optimum  if  the  step  size  is  too  small.   The 
current  step  sizes,  although  generally  too  large,  have  per- 
formed well  for  the  test  craft  above.   The  probable  rationale 
is  that  the  large  step  size  permits  a  search  over  a  wider, 
more  extreme  range  of  system  values,  but  in  the  event  that 
the  step  is  too  large,  the  step  size  will  rapidly  be 
reduced  to  an  acceptable  level. 

The  lack  of  experimental  data  and  analytical  techniques 

has  been  noted  in  several  areas.   The  nacelle  design  is  a 

case  in  point.  Although  nacelles  have  been  studied  exten- 
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sively  by  the  aeronautics  industry,  the  application  of 
nacelles  to  a  water  environment  has  been  neglected  and 
consequently,  design  considerations  have  not  been  defined. 
The  design  of  cascaded  corners  has  somewhat  more  experimental 
data  for  loss  estimation,  but  the  scatter  of  the  data  and 
the  lack  of  definitive  design  information  render  elbow 
design  a  vague  approximation  at  best.   Three  dimensional 
diffusers  are  likewise  an  unexplored  domain.   The  effect  of 
changes  of  ducting  shape  on  losses  is  not  treated  in  the 
literature  for  straight  pipe.   Loss  estimation  for  changing 
shape  in  elbows  is  treated,  but  only  for  rectangular  forms. 
It  is  apparent  that  the  prodram  results  contain  a  great  deal 
of  information,  but  it  is  also  evident  that  the  interpretation 
of  the  results  can  be  improved  by  a  more  thorough  study  of 
runs  for  different  craft,  configurations,  and  missions.   A 
study  of  the  effects  and  interactions  of  such  changes  will 
provide  a  finer  understanding  of  the  import  of  the  results. 
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CHAPTER  4  -  CONCLUSIONS  AND  RECOMMENDATIONS 

A  method  of  optimizing  a  waterjet  propulsion  system 
for  hydrofoils  has  been  developed  and  tested  with  a  four 
hundred  ton  hydrofoil  at  a  cruise  speed  of  forty  two  knots. 
Other  craft  have  verified  the  importance  of  the  selection 
of  the  gas  turbine  for  a  particular  craft.   It  is 
recommended  that  a  version  of  this  program  be  adapted  to 
determine  the  optimum  displacement  for  a  given  engine. 

Although  originally  considered  important,  the  amount 
of  diffusion  in  the  strut  should  be  the  minimum  required  to 
avoid  cavitation  in  the  ducting  or  the  minimum  required  by 
the  pump,  whichever  is  greater.   The  design  of  the  nacelle 
need  not  be  at  the  limits  of  cavitation,  but  should  be 
integrated  into  the  total  system  design  to  provide  a  lesser 
weight  propulsion  system.   It  is  recommended  that  a  study  of 
nacelle  design  be  initiated  to  determine  the  interaction 
of  the  nacelle  and  the  propulsion  system. 

The  initial  values  of  the  system  parameters  may  affect 
the  operation  of  the  optimization  procedure.   A  determination 
of  the  cause  of  this  lack  of  stability  should  be  made.   The 
magnitude  of  the  step  sizes  of  the  parameters  has  little 
effect  on  the  program  provided  they  are  sufficiently  large. 

A  general  lack  of  data  and  design  information  is 

prevalent  for  fluid  flow  in  closed  conduits.   This  is 

especially  noted  for  ram  inlet  design,  cascaded  corners, 

three  dimensional  diff users,  and  changes  of  pipe  shape. 
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The  program  provides  a  unique  facility  for  the 
preliminary  designer  to  rapidly  ascertain  the  optimum 
water jet  propulsion  system  design  and  performance.   Many 
results  should  be  studied  to  deduce  the  complete  real 
world  interpretation  of  a  computer  generated  design. 
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APPENDIX  A  -  OPTIMIZATION 

The  selection  of  an  optimization  procedure  represents 
a  key  consideration  in  the  overall  design  method.   For  one 
thing  it  may  require  computations  to  be  made  solely  for  its 
usage,  as  in  gradient  methods.   Then,  too,  it  may  require 
a  supervisory  program  to  tell  it  what  to  do,  as  in  random 
search  techniques,  or  an  adaptive  direction  capability  may 
be  self-contained  (ref .  18) .   It  is  not  the  intent  of  this 
investigation  to  find  or  invent  the  "optimum"  optimization 
procedure,  but  the  significance  of  its  usage  should  not  be 
minimized. 

The  general  consideration  of  hov;  to  optimize  has  a 

wide  variety  of  tools  to  aid  in  the  achievement  of  the 

answer.   One  of  the  more  common  methods  utilized  is  some 

variant  of  the  gradient  search  (ref.  19).   The  rationale 

is  an  intuitive  one,  based  on  the  notion  that  if  the  goal 

is  the  attainment  of  a  maximum,  for  example,  the  quickest 

way  to  go  is  in  the  direction  of  the  steepest  slope.   This 

is  not  unlike  a  blind  man's  method  of  climbing  a  mountain. 

Like  the  blind  man,  the  gradient  search  can  encounter 

difficulties  by  following  this  philosophy.   The  assumption 

of  a  finite  slope  in  the  area  of  interest  not  only  requires 

the  absence  of  steps  on  the  hypersurface ,  but  imposes  the 

more  stringent  limitation  that  the  slope  be  calculable  and 

containable  within  the  computer.   This  somewhat  limits  the 

types  of  systems  that  can  be  optimized  by  this  technique  and 
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and  exacts  a  stiff  penalty  for  some  numerical  inaccuracies. 
Further,  assuming  that  the  system  to  be  optimized  permits 
slope  calculations,  it  is  not  known  a  priori  that  such 
calculations  provide  sufficient  information  to  compensate 
for  the  time  spent  in  making  them.   In  the  present  system 
design,  discrete  steps  are  virtually  guaranteed  due  to  the 
different  pumps  and  gearboxes  used.   Since  computation  of 
the  propulsion  system  weight  is  not  done  by  means  of  an 
analytic  function,  slope  calculations  become  as  lengthy  as 
the  weight  calculations  themselves.  For  these  reasons, 
gradient  search  techniques  are  rejected  as  a  possible 
optimization  tool  for  this  design  problem  (ref.  20). 

On  the  opposite  end  of  the  optimization  spectrum  lies 
the  random  search.   The  principal  advantages  of  a  random 
search  are  the  ease  with  which  it  m.ay  be  understood,  hence 
used  and  programmed,  and  the  simplicity  of  its  computation. 
No  slopes  or  step  sizes  need  to  be  calculated  or  estimated, 
merely  the  values  of  the  function  at  the  points  of  interest. 
This  is  also  its  major  disadvantage. 

The  random  search  also  has  no  greater  likelihood  of 

finding  a  local  optimum  than  it  does  of  finding  the  global, 

and,  hence  desired  optimum.   Once  the  capability  of 

"directionability  is  gained,  local  optima  are  just  as  likely 

to  be  found  as  the  global  optimum  as  is  the  actual  optimum. 

This  is  true  of  all  directed  searches.   The  only  way  to 

avoid  this  unpleasant  result  is  to  use  random  starting  points 
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or  to  use  a  systematic  search  over  the  entire  hypersurface. 
A  unimodal  surface,  of  course,  does  not  present  this 
problem. 

Since  the  random  search  utilizes  no  other  additional 
information  than  the  function  value,  it  is  highly  in- 
efficient with  the  result  of  greater  time  spent  in  the 
search.   All  optimization  is  done  in  comparison  to  the 
previous  best.   Thus  the  simplicity  of  its  procedure  is 
gained  only  at  the  sacrifice  of  the  large  number  of  evaluations 
that  must  be  made  to  assure  a  reasonably  high  certitude  of 
attaining  the  optimum.   This  becomes  especially  significant 
when  the  time  of  calculation  of  a  single  function  value  is 
lengthy.   In  the  design  and  performance  estimation  of  a 
propulsion  system,  a  lengthy  calculation  is  likely  to  be  the 
case,  in  part  due  to  the  iterative  loops  required  (refs.  21, 
22)  . 

The  solution  is  clearly  somewhere  in  the  middle  of 
these  two  approaches,  utlizing  the  best  features  of  each  and 
minimizing  their  disadvantages.   It  should  have  the  ability 
to  discern  direction  without  the  requirement  of  additional 
calculations.   One  such  method  is  the  pattern  search  utilized 
in  this  design  (ref.  23).   The  basis  of  this  method  is  akin 
to  that  of  the  gradient  search  but  is  modified  in  its  ap- 
proach.  Again  based  on  intuition,  the  pattern  search  depends 
on  the  supposition  that  the  search  effort  should  be  expended 

in  the  direction  of  improvement.   Note  that  this  is  a 
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qualitative  judgment  rather  than  a  quantitative  one,  as  with 
gradient  methods.   The  criterion  for  a  sense  of  direction 
lies  in  the  success  of  an  evaluation,  not  in  its  measure 
of  success,  which  in  itself  would  be  somewhat  arbitrary. 
Thus,  while  gaining  adaptive  information  about  the  location 
of  an  optimum,  the  disadvantages  of  the  gradient  search  are 
neatly  bypassed  and  the  random  search  is  effectively  uti- 
lized in  a  more  efficient  manner. 

Following  the  flow  diagram  of  Figure  8,  the  pattern 
search  uses  the  startina  point  as  the  initial  basepoint. 
From  the  basepoint  exploratory  moves  are  made.   That  is, 
each  coordinate  is  varied  from  the  basepoint  position  the 
amount  of  the  step  size  and  the  function  evaluated  at  that 
point.   Each  exploratory  move  is  examined  for  possible 
improvement  in  the  function  value.   If  improvement  cannot 
be  made  in  any  direction,  the  step  size  is  decremented  and 
the  exploratory  moves  resumed  from  the  basepoint.   If  a 
position  is  located  with  a  better  function  value  than  the 
value  at  the  basepoint,  the  second  type  of  move,  a  pattern 
move,  is  made.   A  pattern  move  is  a  simultaneous  change  in 
all  dimensions  of  length  equal  to  the  distance  from  the 
current  basepoint  to  the  present  position.   If  improvement 
is  made  in  either  the  pattern  move  or  the  sequential  explor- 
atory moves,  another  pattern  move  is  made.   Lack  of  improve- 
ment resets  the  basepoint  at  the  position  of  the  best 

function  value  and  enters  the  cycle  again  with  exploratory 
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moves  from  the  current  basepoint. 

This  ability  to  seek,  find,  and  establish  a  pattern 
gives  this  search  an  ability  to  move  rapidly  in  the  best 
direction.   It  also  has  the  capability  of  shifting  direc- 
tion or  modifying  the  pattern  when  required.   Step  changes 
in  the  hypersurface  present  little  problem  since  slopes  are 
not  utilized.   With  these  features  the  pattern  search  can 
move  quite  readily  along  steep  ridges  or  narrow  valleys 
which  are  common  "Waterloos"  of  other  search  methods, 
especially  gradient  methods. 

Perhaps  it  should  be  mentioned  here  that  the  intent 
is  not  to  totally  degrade  the  gradient  search  which  has  its 
own  attributes  and  an  important  usefulness  in  its  own  right. 
It  is  not  without  shortcomings,  however,  that  might  present 
a  serious  detriment  to  the  user  who  is  looking  for  a 
generalized  search  method. 

The  original  form  of  the  pattern  search  as  proposed  in 

reference  23  has  been  modified  to  conform  to  the  needs  of 

this  design  problem.   An  additional  modification  has  been 

made  that  is  considered  an  improvement  on  the  initial 

scheme.   During  the  exploratory  move  phase  of  the  search, 

each  dimension  was  initially  incremented  by  a  step  size  and 

the  function  evaluated.   If  the  move  was  a  success,  the  point 

and  its  associated  value  were  saved  and  the  next  dimension 

explored.   Lack  of  success  meant  decrementing  the  position 

by  a  step  size  and  evaluating  the  function.   Again  a  success- 
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ful  move  was  followed  by  retaining  the  position  and  value  and 
exploring  the  next  dimension.   An  unsuccessful  move  was 
ignored  and  the  next  dimension  incremented.  The  modification 
implemented  is  to  vary  the  dimension  first  in  the  direction 
of  the  last  improvement  rather  than  simply  incrementing  it. 
Thus  anticipating  success,  time  would  be  saved  by  not 
checking  a  poor  direction.   The  resulting  advantage  of  this 
modification  is,  of  course,  very  dependent  on  the  function 
being  evaluated,  but  is  in  keeping  with  the  philosophy  of 
the  search  and  is  accomplished  at  a  very  small  cost  of  time. 
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APPENDIX    B    -    NACELLE 

The   details   of  this    section  may  be    found  in   reference 
9,    from  which   the    following   comments    are   drawn. 

The   selection   of   the   system  parameter  values   sets   the 
ratio  of  the   inlet  diameter  to  the   maximum  external   diameter, 
DIDM.      This   effectively   limits   that  maximum  amount  of   diffu- 
sion  that  may  be   carried  out   in  the   nacelle.      The  present 
design   assumes   that   the  maximum  amount   of   diffusion    allowed 
in   the  nacelle   duct   is   utilized.      The   nacelle    is    first 
checked  to  ensure    that   it   does  not   cavitate   externally. 
Prediction   of  external   cavitation   rejects   that  selection   of 
system  parameter  values.      If   cavitation    is   not  predicted, 
the   inlet  is   checked   for  cavitation   at  the  point  of  maximum 
flow,   which  generally  occurs   in  the   take-off  condition.      If 
cavitation  was   predicted  internally,    auxiliary   inlets   are 
sized  to  reduce   the   inflow   rate   to   just  below  the  point  of 
cavitation.      The   length  of  the  nacelle    is    then   suboptimized 
to  be   the    length   at  which  the   total  power   loss   due   to  drag 
and  diffusion   in   the   duct  is   minimized.      The   total   head 
loss,   made   up  of  the    lip   loss    and  the   diffuser  losses  may  then 
be   calculated.      Likewise,    the   drag  coefficient  may  be 
calculated   from  Hoerner's   approximation    for  stramlined 
bodies    (ref.    10),    based  on  wetted  surface, 

3 

Cd=Cf    (l+1.5(fn)7+7(Dm)3) 
u        s  In  Ln 
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APPENDIX  C  -  ELBOWS 

All  elbows  in  the  ducting  are  treated  in  identical 
fashion,  although  individually  each  elbow  may  be  very 
different.   The  inlet  and  outlet  areas  for  any  bend  are 
assumed  to  be  the  same.   An  elbow  is  considered  fully  defined 
by  the  knowledge  of  four  quantities:  width,  depth,  crOSS- 
sectional  area,  and  radius  ratio.   The  radius  ratio,  -  , 
is  the  ratio  of  the  radius  of  the  centerline  of  the  bend  to 
the  duct  radius,  as  depicted  in  Figure  9.   Since  the  flow 
rate,  Q,  has  been  previously  determined,  the  first  three 
items  are  sufficient  to  calculate  the  aspect  ratio,  shape, 
and  the  mean  velocity  of  the  fluid  in  the  elbow. 

For  Re (^ ) 2  >  91,  Ito's  expression  for  the  head  loss 
coefficient  in  a  circular  bend  is  (ref .  11) : 

K^.OC^laeRe-0-17^)0-8''  (C.I) 

a  is  a  numerical  factor  deoendent  on  6  and  -  .   For  the 

RO 

case   of   6   =    90°    and  |Q<19.7 

a=0.95+17.2(-    )-1«96  (C.2) 

RO 

The    loss    coefficient  then   is 

Kt=0.0024ie[0.95  +  17.2(go)-1-96]Re""°-  17(f0)°'8,t 

For   a   given   duct   and  considering    G,    V,    and   R  constant, 
partial   differentiation   of  K.    with   respect    to   RO    leads   to 
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9K,.  \       R 


3R0  =    0       ~>    Wopt      *    4'3 


Ci 
Since   K.     =    -~  +   C^RO   for   these    conditions,    this    result 

verifies    the    intuitive    feeling   that   a  minimum  exists    due    to 

a   trade-off  between    friction    and  bend   losses   as   the   number 

of  guide    vanes   in   the   elbow  varies,    as   depicted   in   Figure    10, 

If  the   given   radius    ratio   is   greater   than   one,    the 

losses   in   the   elbow  may  be    reduced  by   the    use    of   splitters 

which   aid  in   the  prevention   of   separation.    In   this   context, 

splitters    are   defined   as   guide   vanes   of   unequal   length, 

unequal    spacing   and  possibly   unequal   curvature    (ref.    12). 

Turning  vanes,    on   the   other  hand,    are    always   equal   in    length 

and  curvature,    although  the    spacing  between   vanes    is   not 

necessarily  equal.      If   the    splitters    are    considered  in   order 

to   subdivide   the   elbow  into  smaller  elbows,    the   minimum 

loss  will   occur  when   the    ratio   of   the    inside    radius   to  the 

outside    radius   of   any   of   the    subdivided  elbows   is   constant. 

Referring   to   Figure    6,    this   may  be   written   as 

(ri,   =    (|i)Vn 

rO    A 

where  n  is  the  number  of  subdivided  elbows,  i.e.  the  number 

of  splitters  +  1.   Figure  11  depicts  a  duct  with  one 

splitter.   Clearly, 

r^  =  R  -  RO 

rQ  =  R  +  RO 
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Hence,  the  desired  radius  ratio,  4.3,  may  be  substituted  for 

(roa+rja)   and  the  resulting  equation  solved  for  n,  or, 
(rOa"ria) 


n= 


Ln 

*--l 
RO 

=    2.12    Ln 

S-+1 

RO    t 

RO 

Ln 

[4.3- 

^] 

*--l 
RO 

L4.3+. 

(C4) 


The  number  of  splitters  determines  the  spacing.  The  head 

loss  coefficient  may  be  calculated  for  each  subdivided 

elbow  and  the  results  averaged  over  the  area  which  will 

yield  the  average  head  loss  coefficient  for  the  elbow. 

For  R  <  1,  equation  C.4  has  no  meaning.  For  this  case 
RO 
the  total  head  loss  coefficient  was  computed  as  the  sum  of 

the  friction  loss  coefficient  and  the  loss  coefficient  from 

Figure  12  which  is  an  empirical  graph  for  thin,  circular  arc 

turning  vanes  from  reference  12. 

In  either  case  the  results  are  valid  up  to  a  Reynolds 

number  below  that  expected  in  this  design.   The  head  loss 

coefficient  including  friction  is  then  corrected  for  the 

higher  Reynolds  number  by  (ref.  13) 


K2  "  2k  K! 
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where   by   curve    fit, 

X  =    1.0057-0. 16892Ln(Re   x   10"5)    +0 . 0145885 [Ln (Re   x    10~5)]2 

-0.0016948[Ln(Re    x   10"5)]3  (C.5) 


-55- 


APPENDIX  D  -  STRUT 
The  head  losses  in  the  strut  diffuser  are  caused  by 
friction  and  the  area  expansion  which  thickens  the  boundary 
layer.   The  friction  coefficient  is  computed  in  the  same 
manner  as  for  straight  pipe.   The  expansion  loss  coefficient 
is 

Ke  =  C(l-  |i)2  (D.l) 

C  is  dependent  on  the  diffuser  equivalent  angle,  20,  and 

is  interpolated  from  the  graph  in  Figure  13,  which  is  for 

a  two- dimensional,  straight-walled  diffuser.   For  the  flow 

area,  the  thickness  to  chord  ratio  of  the  strut  is  held 

constant  while  the  area  is  increased.   Hence  the  assumption 

of  two  dimensional  diffusion  is  not  strictly  valid  but  is 

considered  a  conservative  estimate. 

The  strut  dissipates  energy  externally  by  a  drag  term 

attributable  to  the  body  and  a  term  resulting  from  the  spray. 

The  drag  coefficient  from  Hoerner  (ref.  10)  is  based  on  the 

planform  area  and  is  approximated  by 

C.   =  2Cf  [l+2(t/c)+60(t/c)"]  (D.2) 

us      s 

The  spray  drag  coefficient,  likewise  based  on  planform  area 
is 

C-,   =  .03(t/c)  (D.3) 

SP 
This  approximation  is  valid  for  Froude  number,  based  on  the 

chord  length,  greater  than  three  (ref.  14) . 
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APPENDIX  E 

Pump  Inlet  Piping 

The  ducting  from  the  hull  elbow  to  the  pump  inlet  is 

a  straight  pipe,  a  9  0°  elbow  and  a  transition  piece.  The 

length  of  the  straight  pipe  depends  on  the  number  of  prime 

movers.   In  equation  form, 

Length  =   Beam 

Number  of  prime  movers   +    1 

The   transition  piece    is   shaped   according   to  the   number  of 
pumps   under   consideration.      Since   there    are   two   struts 
ducting  water  into  the   hull,    assuming   a   canard   foil 
arrangement,    and  only   one,    two,    or   four  pumps   are    consi- 
dered,   the   transition  piece  will  be   a   junction,    straight 
pipe,    or  diverging  branch. 

All   straight  pipe    loss    calculations    are   made   using   the 
standard  Darcy-Weisbach  equation 

h  =    f L  V^  (E.l) 

De2g 

The  Moody   friction   factor,    f,    is   computed  using  the   closed 

form  approximation   of  Techo  et   al    (ref .    15) , 

f  =    [.86859    Ln(Re/1.964Ln(Re)    -    3.8215))]"2  (E.2) 

R     is   based  on   the  mean  hydraulic  diameter  which  is   defined 

as 

De   =    4A  (E.3) 

P 

This  approximation  for  the  friction  factor  is  accurate  to 

within  .1%  over  the  range  of  Reynolds  number  from  10   to 
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2.5   x    108. 

For  both  the    junction   and  divergence,    the   average 
velocity   is   assumed  to  be   the    same   before    and   after  the 
transition  piece    and  the  piping  is   symmetrical  about   the 
single   pipe    as   shown   in   Figure    14.      This    requires    that  the 
smaller  piping   cross-sectional   area   is   1/2   that  of   the 
larger  piping. 

The  mixing   loss    coefficient  in    the    junction,    without 
friction    loss,    is   computed   from 

Kj   =    1   +    X    -    2cos(f  ($)  )  (E.4) 

where    f(<f>)    =    1.4$   -    .00583(J)2 

and   <f>   is    the    angle    the    small  branch  makes  with  the    larger 
duct    (ref.    16).      In   the    cases   examined   for   this   design,    <j> 
was   set  to   zero   resulting  in   only   friction   losses   in   the 
junction. 

The   divergence    loss    factor  X  was   computed   from  Figure 
15  which  is   adapted   from   reference    16.      The   angle  of   the 
divergence,    of  course,    depends   on   the  pump   length.      It 
should  be  noted   that   the   divergence   configuration   requires 
another  elbow,    equal  in   angle   to  that  of  the   divergence, 
to  enter  the  pump. 

In   all  three   transition  pieces   the   size  was    set  as 

required,    depending  on   the  pump   length.      If  the  pump's 

length   requires   a   zero  or  negative   transition  piece 

length,    the   transition  piece   is  neglected  and  the   losses 

computed  only    for  the    athwart ships    length   and  the  pump  elbow 
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which   turns   the   ducting   aft. 
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APPENDIX    F 

Pump 

The   details   of   this   appendix  may  be    found  in   reference 

17,    from  which   the    following   comments    are    drawn. 

There    are    two  basic   types   of  pumps   considered   in    this 

design:      a  multistage,    double    impeller   centrifugal   design 

and  an   axial  pump  with   an   inducer  impeller.      The   variations 

occur  in    the   number  of  stages   used.      The  pump   design 

utilizes    the    flow   rate    and  Thoma's    criterion   is    defined  as 

a     =   NPSH  (F.l) 

Pump  Head 

where  the  net  positive  suction  head, 

NPSH  =  v£  +  ha  -  hv  -helev  -  Ahi  (F.2) 

2a 

and 

Pumo  Head  =  V?   -  V2     +  hQlo„  +   Alv  (F.3) 

— j  o  eiev  t 

2g 

The   axial  pump's    first   stage   is   an   inducer  impeller 
with   up   to   two   additional   stages   available    if   the    amount 
of  head   is   sufficiently   high.      Assumption   of   a   specific 
speed  yields    the    design  pump   speed   and   approximation 
to   the   off-design  pump    speed,    which   in    this    case    is   the 
cruise    condition    since   the   pump   is  basically  designed 
for  take-off   conditions.      The   RPM  of   the  prime   mover  deter- 
mines  the   gear  ratio   and  calculation   of   the    shaft  horse- 
power permits   an  estimation  of   the   gear  weight.      Finally, 

a   fuel  weight  may  be    calculated   for   the    range    specified   in 
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the   cruise   condition. 

The    centrifugal  pump   is    designed  in    the    same   way  except 
that   the   number  of  double   suction   impellers   is   increased, 
thus   increasing   the   RPM,    until   the    combined  weights   of  the 
pump,    dry   and  wet,    gearbox  and   fuel  are   at   a  minimum. 

The    two   resulting  pump   designs   are    then   compared  on 
the   total  pump,    fuel,    and  gearbox  weights    and   the    least 
weight   system  is    selected   for  the   propulsion  plant. 
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APPENDIX   G 
Nozzle 
For  the  purpose   of  calculation   the  nozzle  was   consid- 
ered a   straight-walled  nozzle.      Assuming   the    losses    are  of 
the  .form 

h     =    f  L  V^  (G.l) 

n  D   2g 

and  noting  the  friction  factor  is  approximately  constant 

and  the  velocity  essentially  uniform,  this  may  be  written 

Ah     =    f  AL  V2 
n  D  2g 

or 

L-i  2 

hn   =    f      P    1   \T    dL 

0      D   2g 
Substituting        D  =   D.    +  Ltan(i|0 


V  =    4Q 
results   in 


ttD2 


16  V 


hn  =    f        J/2    [1-(1+Dt)*]    12  (G-2) 

Dt"Dj  crj  2? 

which  is   used   for  the   head  loss  estimation   in   the   nozzle 

when   D        >    D.      >    1,     (ref.    2).      When   D.      <    1,    the   head  loss 

t     j  j 

is  assumed  to  be  1.5%  of  the  dynamic  head. 

In  the  event  that  the  distance  of  the  pump  exit  from 

the  transom  is  greater  than  2  times  the  height  of  the  pump 

centerline  above  the  keel,  the  nozzle  is  assumed  to  exit 

through  the  underside  of  the  craft.   The  latter  distance 
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is  that  if  the  nozzle  were  directed  along  an  axis  at  45° 
from  the  vertical. 

In  every  instance,  however,  the  nozzle  is  assumed  to 
be  depressed  at  the  optimum  angle  to  reduce  the  lift 
required  which  would  allow  the  use  of  smaller  foils  and 
hence  reduce  drag.   As  before, 

T  =  pQ(V.cos(B)  -  VQ)  (1) 

Rewritten   this   is   also 

T   =    (W   -   T    )    D  (G.3) 

L 

where 

Tv  =    pQV.sin(S) 

Equating  equations  (1)  and  (G.3),  and  solving  for  V.  gives 

V.  =  (PQVQ  +  WD) [cos  ($)  +  Dsin  (B)]"1 
3  L  L 

The  minimum  energy  input  to  the  jet  occurs  when  V.  is  at  a 

minimum,  or, 

3V-:  =  0  =>  (pQV0  +  WD)(Dcos($)  -  sin(3)) 

9BJ        L  L  =  0 

(cos (6)  +  Dsin(B) P 
L 

=>  3  =  TAN_1r  D  i 

L  (G.4) 

which  is  the  optimum  angle  for  nozzle  depression  (refs.  4,6) 
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FIGURE    9  .      ELBOW 


-65- 


K, 


R/4.3  RO 


FIGURE  10    MINIMUM  ELBOW  LOSS  COEFFICIENT 


-66- 


V 
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RETURN  TC  CALLING 
PROGRAM 


Label 


PARMS 


DRAG 


FLOW 


Variable 

Name 


VJVO 
VIVO 
DIDM 


LIST  OF  SYMBOLS 
USED  IN  COMMON 
Subscript 


TDRAG(I) 

I 

1 

2 

3-5 

STRTD(I) 

POD(I) 

SPRAY  (I) 

REST  (I) 

VO(I) 

TRIM  (I) 

Q(D 


Description/Usage 

System  parameters 

Jet  velocity  ratio 

Inlet  velocity  ratio 

Inlet  diameter  to 
maximum  nacelle  diameter 
ratio 

Drag  estimates 

Total  craft  drag  at 
condition  I,  pound  force 

Cruise  condition 

Take-off   condition 

Evaluation   conditions 

Strut  drag  at   condition 
I ,   pound   force 

Nacelle   drag   at   condition 
I ,    pound   force 

Spray   drag   at   condition 
I,    pound   force 

TDRAG  (I)  -SPRAY  (I)  -STRTD  (I) 
-POD(I),   pound   force 

Craft   velocity  at   condition 
I,    foot  per   second 

Craft  trim  angle  at  con- 
dition  I,    degree 

Flow  characteristics 

Volume  flow  rate  at  con- 
dition I,  cubic  foot  per 
second 
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Label 


ELBW 


CHARS 


Variable 

Name 

AIN 

AJET 

AREA(ICOMP) 

VJ(I) 
VI  (I) 


Subscript 


XK(IELB) 


IELB 

1 
2 
3 
4 


RO(IELB) 


THATA(IELB) 


WIDTH 


DEPTH 


TYPE (3, I TYPE) 


WGTS (LS,LC) 


Description /Us age 

Inlet   area,    square    foot 

Jet    area,    square    foot 

Ducting   cross-sectional 
area  at   station   ICOMP, 
square    foot 

Jet  velocity   at  condition 
I,    foot  per  second 

Inlet  velocity   at  con- 
dition  I,    foot  per   second 

Elbow  data 

Radius  ratio  of  elbow 
IELB 

Strut  elbow 

Hull  elbow 

Pump  elbow 

Divergence  elbow 

Duct  radius  at  elbow 
IELB,  foot 

Angle  of  bend  of  elbow 
IELB,  degree 

Width  of  duct  at  elbow 
inlet 

Depth  of  duct  at  elbow 
inlet 

Contains  name  of  elbow 
shape  ITYPE 

System  characteristics 

Weight  of  component  LC, 
pound 
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Variable 

Label 

Name 

Subscript 
LS 

Description/Usage 

1 

Structure 

2 

Water  (or  fuel) 

LC 

1 

Nacelle 

2 

Strut  elbow 

3 

Strut  diffuser 

4 

Hull  elbow 

5 

Athwart ships  length 

6 

Pump  elbow 

7 

Transition  piece 

8 

Pump 

9 

Nozzle 

10 

Reduction  gear 

11 

Fuel 

12 

Prime  mover 

13 

Lift  from  nozzle 
depression 

14 

Total  weight 

15 

Soare  location 

CGS(LG,LC) 

Centers  of  gravity 

LG 


component  LC,  excluding 
fuel,  gearbox,  or  prime 
mover,  foot 


1  Structure  vertical  center 

gravity  (from  keel) 
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Label 


Variable 

Name 


Subscript 
2 


DELH(I,ICOMP) 


ICOMP 

8 

9 


10-15 


CGSX 


CGSZ 


H20 


TEMP 

PV 

RHOW 


Description/Usage 

Structure    longitudinal 
center  of   gravity 
(from   transom) 

Water  vertical   center  of 
gravity,     (from  keel) 

Water  longitudinal  center 
of   gravity,     (from  transom) 

Total  head  loss  up  to  and 
including  component  ICOMP 
at   condition   I,    foot 

Same    as   defined  except 

Nozzle   head   loss   only 

Total   system  head   loss 
(including  elevation) , 
DELH ( I , 9 ) =DELH ( I , 8 ) + 
DELH    (1,7) 

Spare    locations 

Total    longitudinal   center 
of  gravity  of  system 
excluding  gearbox,    prime 
mover  and  fuel,    from 
transom,    foot 

Total   vertical   center  of 
gravity   of   system  exclud- 
ing gearbox,   prime 
mover  and   fuel,    from  keel, 
foot 

Sea  water  (3.5%  salinity) 
properties 

Temperature,    degree 
Fahrenheit 

Vapor  head,    foot 

Density,    pound   force 
second   squared  per   foot 
to  the   fourth 
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Label 


TOLER 


PSUB 


Variable 

Name 

GNU 


HA 

DELTA 

GERAT(NSTG) 
SHP(I,NSTG) 

RPM(I,NSTG). 

PERF(L,IENGN) 


Subscript        De  s  crip ti  on/Us  age 

Viscosity,    foot 
squared  per   second 

Atmospheric  head,    foot 

Check    for  optimum  system 

Pump   system 

Gear   ratio   required   for 
pump  NSTG 

Shaft   horsepov/er   re- 
quired  at   condition   I    for 
pump   NSTG,    horsepower 

Axial    speed  of  pump   NSTG 
at   condition   I,    revo- 
lution per  minute 

Prime  mover  IENGN 
Characteristics    (L) 


L 

1 


4 
5 


ETAP(I,NSTG) 


Maximum  normal  horsepower 
at  design  speed,  horse- 
power 

Maximum  intermittent 
horsepower  at  design 
speed,  horsepower 

Specific  fuel  consump- 
tion (SFC)  at  design  speed 
and  maximum  normal  horse- 
power, pound  fuel  per 
horsepower  hour 

Design  speed,  RPM 

Prime  mover  weight,  with- 
out auxiliaries,  pound 

Efficiency  of  pump  NSTG 
at  condition  I 
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Label 


Variable 

Name 


SHIP 


NACELLE 


DISP 

RANGE 

BEAM 

HS 

HE 

HCL 

XLS 

XLPE 
XLP 

DRAT 
DM 

AI 

AIAUX 

ELEXT 
ELENT 


Subscript 

I 

7 

6,8 


Description/Usage 

Same  as  defined  except 

Reduction  gear 

Spare  locations 

Craft  characteristics 

Craft  displacement,  pound 

Endurance,  nautical  mile 

Beam,  foot 

Depth  of  submergence  of 
foil,  foot 

Height  of  pump  center line 
above  keel,  foot 

Height  of  pump  centerline 
above  mean  water,  foot 

Distance  of  centerline 
of  strut  root  from  tran- 
som, foot 

Distance  of  pump  exit 
from  transom,  foot 

Length  of  pump,  foot 

Nacelle  characteristics 

Diameter  ratio,  DI/DM 

Maximum  external  diameter, 
foot 

Inlet  area  per  nacelle, 
square  foot 

Auxiliary  inlet  area  per 
nacelle,  square  foot 

Length  of  forebody,  foot 

Length  of  lip,  foot 
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Label 


CONST 


STRTC 


INDEX 


Variable 

Name 

ELAUX 


ELDIF 
ELN 


PI 
G 

RHOD 


TC 

T 

C 

Tl 

CI 

CFM 


IEVAL 


Subscript   Description /Usage 

Length  due  to  auxiliary 
inlet,  foot 

Length  of  diffuser,  foot 

Length  of  nacelle,  foot 

Constants 

3.14159265 

Acceleration  of  gravity, 
32. 174, foot  per  second 
squared 

Density  of  steel,  480, 
pound  per  cubic  foot 

Strut  characteristics 

Thickness  to  chord  ratio 

Thickness  at  root,  foot 

Chord  at  root,  foot 

.  Thickness  at  tip,  foot 

Chord  at  tip ,  foot 

Chord  at  flying  water- 
line,  foot 

Indices  for  program 
control 


<0 


=0 


>0 


Design/evaluation  index 

Evaluate  at  | IEVAL |  points, 
no  design 

Design  at  cruise  and  take- 
off conditions  only 

Design  at  cruise  and  take- 
off conditions  and  evaluate 
at  | IEVAL |  points 
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Label 


Variable 

Name 

IEQPT 


Subscript   Description/Usage 


ISTRT 


NUMB 


IENGN 


=  0 


>0 


1 
3 


2+|lEVAL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 


12 


No  equipments  or  config- 
uration entered;  total 
system  design/evaluation 

Specific  equipments  or 
configuration  entered 
(not  yet  implemented  into 
program) 

Initial  program  condition 

Design  and/or  evaluate 

Evaluate  only 

Final  program  condition 

Design  only 

Design   and/or  evaluate 

Prime   mover  type 

TF   35 

TF   40 

Proteus,    1500    rpm 

Proteus,    1000    rpm 

Tyne    1A 

Tyne    1C 

FT12A 

LM    1500 

LM   2500 

FT  4 A- 2 C 

FT4A-12 

FT4C-2 
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Label 


Variable 

Name 

I  TYPE 


ICOMP 


NPUMP 

NGT 

ITABL 

CDRAG 

CSTRT(I) 

CPOD(I) 

CSPRY(I) 

Subscript 

1 
2 

3 
4 

1 
2 
3 
4 
5 
6 
7 
8 
9 


WEGT 


Description/Usage 

Elbow  shape 

Ellipse 

Circle 

Rectangle 

Square 

Component  index 

Nacelle  outlet 

Strut  elbow 

Strut  diffuser  exit 

Hull  elbow 

Athwartships  length 

Pump  elbow 

Transition  piece 

Pump   inlet 

Nozzle  throat 

Number  of  pumps 

Number  of  gas  turbines 

Interpolation  parameters 

Computed  drags 

Computed  strut  drag  at 
condition  I,  pound  force 

Computed  nacelle  drag  at 
condition  I,  pound  force 

Computed  spray  drag  at 
condition  I,  pound  force 

Pump,  gear  and  fuel  weights 


-81- 


Label 


PUMM 


Variable 

Name 

XWD(NSTG) 


XWW(NSTG) 
XWG(NSTG) 
XWF(NSTG) 


QQ(D 

DIS(NSTG) 

D2S (NSTG) 

XNS(NSTG) 

SM(NSTG) 

PLP(NSTG) 
NSTG 


Subscript   Description/Usage 

Pump  NSTG  dry  weight, 
pound 

Pump  NSTG  wet  weight, 
pound 

Gearbox  weight  associated 
with  pump  NSTG,  pound 

Fuel  weight  associated 
with  pump  NSTG,  pound 

Pump  characteristics 

Flow  rate  per  pump  at 
condition  I,  cubic  foot 
per  second 

•  Inlet  tip  diameter  of 
pump  NSTG,  foot 

Exit  tip  diameter  of  pump 
NSTG,  foot 

Specific  speed  of  pump 
NSTG,  cfs  units 

Suction  specific  speed 
of  pump  NSTG,  cfs  units 

Length  of  pump  NSTG,  foot 

Indicator  of  pump  type 


1 
2 
3 

4 


Axial  pump  with  single 
inducer  impeller 

Axial  pump  with  inducer 
and  one  impeller  stage 

Axial  pump  with  inducer 
and  two  impeller  stages 

Centrifugal  pump  with  max- 
imum of  ten  parallel 
double  suction  impellers 

Spare  location 
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Label 


HEAD 


***** 


Variable 

Name 

SHI (NSTG) 


XIM 


HPP  (I) 

HSV(I) 

THOM(I) 

PHI (NSTG) 

WF 

WG 


IFUEL 


Subscript        Description/Usage 

Head  coefficient   of  pump 
NSTG 

Number  of  parallel   double 
suction   impellers    for 
centrifugal  pump 

Number  of  impellers,    not 
including   inducer,    for 
centrifugal  pump 

Pump    flow   characteristics 

Pump   head   at   condition   I, 
foot 

Net  positive    suction   head 
at  condition   I,    foot 

Thoma's    cavitation   index 
at  condition   I 

Flow  coefficient  of  pump 
NSTG 

Working  variable    for   fuel 
weight,    pound 

Working  variable    for  gear 
weight,    pound 

Unlabeled  common   variables 

Logical   variable    for   fuel 
calculation 

.TRUE.  Make    fuel   calculation 

.FALSE.  Do  not  make    fuel   calcula- 

tion 
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Variable 

P ARM ( 3 ) 

ENGN(3,IENGN) 

EHP(I) 

DEL (3) 

DELMIN(3) 

VK(I) 

IPRNT 

DISPL 

IN 

MAX 

XJ 
IK 
IL 
KM 

Kl 

WEIGT 

IM 

DISPL 

I 

IJ 
J 


VARIABLES  USED  IN  H20JT 
NOT  IN  COMMON 

Description 

Working  system  parameters 

Contains  name  of  engine  IENGN  for  output 

Effective  horsepower  at  condition  I,  horsepower 

Initial  step  sizes  of  system  parameters 

Minimum  step  sizes  of  system  parameters 

Craft  speed  at  condition  I,  knot 

Print  data  set  reference  number 

Craft  displacement,  long  ton 

Estimate  of  minimum  number  of  gas  turbines 
required  to  power  craft 

Estimate  of  maximum  number  of  gas  turbines 
needed  to  power  craft 

Working  variable  for  number  of  gas  turbines 

Working  variable  for  craft  condition  I 

Working  variable  for  IN 

Minimum  number  of  pumps  to  be  used  for  XJ 
gas  turbines 

Maximum  number  of  pumps  to  be  used  for  XJ 
gas  turbines 

Propulsion  system  weight,  long  ton 

Working  variable  for  prime  mover  character- 
istic L 

Craft  displacement,  long  ton 

Index  for  do  loops 

Working  variable  for  check  on  power  adequacy 

Index  for  do  loop 
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Variable 
I  SAVE 

N 

XMIN 

MK 


Description 

Working  variable  for  MAX 

Working  variable  for  output  of  prime  mover 
tyoe 

Factor  for  minimum  acceptable  prime  mover 
SHP  operation 

Working  variable  for  prime  mover  character- 
istic L 
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LIST  OF  VARIABLES  USED  IN  FCT 
NOT  IN  COMMON 

Variable  Description 

APOD  (I)        Difference  between  present  and  previous 

nacelle  drag  calculations  at  condition  I, 
pound  force 

ASTRT(I)       Difference  between  present  and  previous  strut 

drag  calculations  at  condition  I,  pound 
force 

ASPRY(I)       Difference  between  present  and  previous  spray 

drag  calculations  at  condition  I,  pound 
force 

C(I)  Margin  factor  for  thrust  at  condition  I 

TSUM  Workina  variable  for  propulsion  system 

weight,  pound 

DSIG  Working  variable  for  SIGMA 

TOTAL  Total  head  less  vapor  head,  foot 

SIGMA         Head  above  vapor  head,  foot 

SUMZ  Working  variable  for  moment  in  vertical  plane, 

due  to  weights,  foot  pound  force 

SUMX  Working  variable  for  moment  in  longitudinal 

plane  due  to  weights,  foot  pound  force 

SUM  Working  variable  for  system  weights,  pound 

ZSUM  Total  moment  in  vertical  plane  due  to  weights, 

foot  pound  force 

XSUM  Total  moment  in  longitudinal  plane  due  to 

weights,  foot  pound  force 

SPD  Working  variable  for  jet  velocity  calcula- 

tion, foot  per  second 

J  Working  variable  for  NUMB 

JDRAG  Indicator  of  acceptable  drag  accuracy 

1     New  drag  calculation  is  within  5%  of  previous 
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Variable 


Description 


ANGLE 
COEF 


KOUNT 

I 

PARM(K) 


DDRAG ( I ) 
DSPRY(I) 
DSTRT(I) 
DPOD(I) 


0 
1 


K 
1 
2 


calculation 

New  drag  calculation  is  greater  than  5%  of 
previous  calculation,  redesign  strut  and 
nacelle 

Optimum  nozzle  depression  angle,  radian 

Factor  for  nacelle  drag  calculation 

Strut  must  be  resized  due  to  cavitation 

Use  nacelle  design 

Working  variable  for  craft  condition  I 

Craft  condition 

Working  variable  for  system  parameters 

VJVO 

VIVO 

Previous  acceptable  total  drag  calculation  at 
condition  I,  pound  force 

Previous  acceptable  spray  drag  calculation  at 
condition  I,  pound  force 

Previous  acceptable  strut  drag  calculation  at 
condition  I,  pound  force 

Previous  acceptable  nacelle  drag  calculation 
at  condition  I,  pound  force 
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LIST    OF   VARIABLES    USED    IN   NACEL 
NOT    IN    CO?lMON 

Variable  Description 

ZK  Percentage   of   auxiliary   inlet  area  permitting 

flow 

SPO  Static  pressure  at  inlet,  pound  per  square 

foot 

PVP  Vapor  pressure,    pound  per   square    foot 

SIGTV  Incipient   turning  vane    cavitation  number, 

referenced   to   diffuser  exit  velocity 

JNUMB  Working   variable    for  NUMB 

I  Craft  condition,    also   index   for  do    loops 

CPEX  Peak   external  pressure    coefficient 

K  Working  index   for  do    loop    for  ELEXT/DM 

values 

J  Working  index  for  do  loop  for  VIVO  values 

XD  Length  of  forebody  to  maximum  diameter  ratio, 

ELEXT/DM 

DIDMX  Maximum  permissible   diameter   ratio 

QI  Flow   rate  per  strut,    cubic   foot  per   second 

DI  Nacelle    inlet  diameter,    foot 

CPIN  Peak   internal  pressure    coefficient 

QIN  Take-off    flow   rate   per  nacelle    through   inlet, 

cubic    foot  per  second 

QC  Take-off    flow   rate   per  nacelle,    total,    cubic 

foot  per  second 

QAUX  Auxiliary    flow   rate    required   to   avoid 

cavitation,    cubic   foot  per  second 

KDEX  Counter   for   iterations    on   diffuser   length 

VR  Velocity   ratio   at  take-off  condition,    based 

on    flow   through   inlet 
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Variable 


Description 


VI 2 

PRL2 

PTI 

SPI 

VIAUX 

PRAUX 
PTAUX 

DYP 

PC 

QDIF 

PHI 

PHS 

X 

D2 

Dl 

ELMAX 

ELMIN 

II 

EL 

DEL 


Inlet  velocity   at  take-off   condition, 
based  on   flow   through   inlet 

Pressure    recovery   coefficient  of   lip   at 
take-off 

Inlet   stagnation  pressure    immediately   aft 
of   lip   at  take-off,    pound  per   square    foot 

Inlet   static  pressure    immediately   aft  of   lip 
at  take-off,    pound  per   square    foot 

Inlet  velocity   in   auxiliary   inlet,    foot  per 
second 

Pressure   recovery  of  auxiliary   inlet 

Stagnation  pressure    inside    auxiliary   inlet, 
pound  per  square    foot 

Net  dynamic  pressure    immediately   aft  of   lip, 
pound  per  square    foot 

Average    inlet   stagnation  pressure   of   combined 
flow,   pound  per  square    foot 

Working  variable    for  diffuser   flow   rate, 
cubic   foot  per  second 

Equivalent   angle    of   forebody,    radian 

Sine   of  equivalent   angle   of    forebody 

Length   of   auxiliary   inlet,    foot 

Diffuser  exit  diameter,    foot 

Inlet  diameter,    foot 

Maximum  permissible    length  of  diffuser,    foot 

Minimum  permissible    length   of   diffuser,    foot 

Working  variable    for  craft  condition   I 

Working  variable    for  diffuser   length,    foot 

Working  variable    for   change    in   diffuser   length, 
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Variable 


Description 


ELD 
ELL 
ELF  AC 

DDM 
XKT 
REL 
RED 
DL 

CDRG 
ANGL 
CDIF 
POW 

POWI 

EM 
AEXN 
REND 
DDIF 

PLOSS 

VAOUT 
SQUAR 
VCRIT 


Working  variable  for  nacelle  length,  foot 

Working  variable  for  nacelle  length,  foot 

Working  variable  for  excess  in  nacelle  length 
due  to  diffuser,  foot 

Average  diffuser  diameter,  foot 

Form  loss  coefficient  of  diffuser 

Reynolds  number,  based  on  nacelle  length 

Reynolds  number,  based  on  inlet  diameter 

Ratio  of  maximum  external  diameter  to  nacelle 
length 

Computed  drag  coefficient 

Equivalent  half  angle  of  diffuser,  degree 

Diffuser  expansion  factor 

Power  loss  due  to  drag  and  duct  loss  of 
diffuser,  horsepower 

Previous  power  loss  calculation  for  diffuser 
length,  horsepower 

Factor  in  wetted  surface  calculation 

Wetted  surface  area,  square  foot 

Reynolds  number,  based  on  inlet  diameter 

Total  pressure  loss  in  diffuser,  pound  per 
square  foot 

Total  pressure  loss  in  nacelle,  pound  per 
square  foot 

Average  exit  velocity,  foot  per  second 

Factor  in  critical  velocity  calculation 

Critical  velocity  in  strut  elbow,  foot  per 
second 

-90- 


Variable 


Descriotion 


VMAX 


RENL 

NL(2) 

ML(2) 

KL(2) 

JL(2) 

IL(2) 

VRT(6) 

XDT(IO) 

PRLT(6) 

SIGI  (I) 

PTO(I) 

XDTT(J) 

CDUMX (K) 

CDUMY(L) 

CD  (I) 

DIDMT(J) 

VRTT(J) 

DLIP(I) 

QO(I) 

VELR(I) 


Maximum  velocity  at  nacelle  exit,  foot  per 
second 

Reynolds  number,  based  on  nacelle  length 

Interpolation  parameters 

Interpolation  parameters 

Interpolation  parameters 

Interpolation  parameters 

Interpolation  parameters 

Data  array  of  velocity  ratios 

Tabulated  forebody  length  to  inlet  diameter 
ratios 

Tabulated  lip  pressure  recovery   coefficients 

Free  stream  cavitation  index  at  condition  I 

Stagnation  pressure  at  craft  condition  I, 
pound  per  square  foot 

Tabulated  ELEXT/DM  ratios  for  trim  angle  J 

Dummy  array  of  peak  pressure  coefficients  for 
VIVO  K 

Dummy  array  of  velocity  ratios  for  ELEXT/DM 
L 

Drag  coefficient  at  craft  condition  I 

Tabulated  DI/DM  values  for  ELEXT/DM  J 

Dummy  array  of  velocity  ratios  for  angle  of 
attack  (J) ,  internal 

Lip  loss  coefficient  at  craft  condition  I 

Free  stream  dynamic  pressure  at  craft  condition 
I,  pound  per  square  foot 

Inlet  velocitv  ratio  at  craft  condition  I 
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Variable 
VRTEX(I) 
VRMAX(I) 


Description 

Same  as  VRTT(I)  but  external 

Maximum  permissible  velocity  ratio  at  craft 
condition  I 


FRICT 


Function  statement,  calculates  Moody  friction 
factor  for  smooth  pipe,  based  on  Reynolds 
number 
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Variable 
SHAPE(3,ITYPE) 
THETA(J) 
XLOSS(J) 

ROA(IO) 

RE 

FRICT 

CORR 

REM  AX 

KOUNT 
IELB 

1 
2 
3 
4 

FACTR 

RIN 

ROUT 

RATIO 

XN 


VARIABLES  USED  IN  ELBOW 
NOT  IN  COMMON 

Description 

Contains  name  of  shape  ITYPE 

Data  array  of  elbow  angles 

Data  array  of  elbow  loss  coefficients  with 
thin,  circular  arc  turning  vanes 

Outside  radius  of  splitters,  foot 

Reynolds  number  of  duct 

Function  statement,  calculates  Moody  friction 
factor  for  smooth  pipe,  based  on  Reynolds 
number 

Function  statement,  calculates  head  loss 
correction  factor  for  different  Reynolds 
numbers 

Maximum  Reynolds  number  permitted  for 
splitter  loss  equation 

Working  variable  for  craft  condition  I 

Index  indicating  which  elbow  is  being 
designed/evaluated 

Strut 

Hull 

Pump 

Divergence 

Factor  used  in  splitter  loss  calculation 

Inside  radius  of  bend,  foot 

Outside  radius  of  bend,  foot 

Desired  radius  ratio,  4.3 

Number  of  subdivided  elbows  required  to 
achieve  RATIO=4.3 
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Variable 


Description 


N 
Nl 

SUM 

RIA 
V 

HOT 
AA 

RAD 
XCORR 

XKT 
DIAM 
VOLV 
VOL 

AREA1 
IJ 
IK 
I 

RATEO 

HEAD 


Number  of   subdivided  elbows    used 

Number  of   splitters   corresponding   to  N 

Working   variable   of   sum  of  head  times 
subdivided  elbow   area,    foot   cubed 

Inside    radius   of   subdivided  elbow,    foot 

Average   velocity,    foot  per   second 

Height   of   subdivided  elbow,    foot 

Equivalent   cross-sectional   area  of   subdivided 
elbow  based  on  HGT,    square    foot 

Equivalent   radius   of    subdivided  elbow,    foot 

Ratio  of  head   loss    correction    factors    for 
differing   Reynolds   numbers 

Head   loss    coefficient   of   subdivided  elbow 

Equivalent  diameter,    foot 

Volume   of   splitters,    cubic   foot 

Volume   of   splitters   and  elbow  structure, 
cubic   foot 

Duct   area  per  elbow,    square    foot 

Working  variable    for  shape   determination 

Working  variable    for  shape   determination 

Working   variable    for  number  of   subdivided 
elbows   N 

Ratio  of   inside    radius   to  outside    radius 
of   subdivided  elbows 

Average   head   loss   in  elbow,    foot 


-94- 


LIST  OF  VARIABLES  USED  IN  STRUT 
NOT  IN  COMMON 

Variable  Description 

THET2  Array  of  data  of  equivalent  angle  of  diffuser, 

degree 

EXPAN  Array  of  data  of  expansion  coefficient 

ARATO  Area  ratio  of  strut  diffuser 

VL  Local  maximum  velocity,  external,  foot  per 

second 

SIGMA  Local  cavitation  number,  external 

CM  Mean  chord,  foot 

DEIN  Strut  inlet  equivalent  diameter,  foot 

WIDE  Width  of  duct  at  strut  exit,  foot 

DEOUT  Strut  exit  equivalent  diameter,  foot 

DEAVE  Average  equivalent  diameter,  foot 

STRT  Vertical  strut  length,  foot 

XLONG  Actual  strut  length,  foot 

STAN  Arctangent  of  equivalent  diffuser  angle 

THETA  Equivalent  angle  of  diffuser,  26,  degree 

ECOEF  Diffuser  expansion  factor 

FORML  Diffuser  expansion  loss  coefficient 

VELIN  Average  inlet  velocity,  foot  per  second 

VLOUT  Average  exit  velocity,  foot  per  second 

RES  Strut  Reynolds  number,  based  on  mean  chord 

CDS  Strut  drag  coefficient 

CDSP  Spray  drag  coefficient 

RE  Duct  Reynolds  number,  based  on  inlet  velocity 
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Variable 

PIPEL 

TOTAL 

HEAD 

KOUNT 

HGT 

CGWS 


Description 
Duct  friction  loss  coefficient 
Total  loss  coefficient 
Head  loss  of  diffuser,  foot 
Working  variable  for  craft  condition  I 
Elevation  of  strut,  foot 
Vertical  center  of  gravity  of  duct,  foot 
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LIST    OF   VARIABLES    USED    IN    JUNCT 
NOT    IN    COMMON 

Variable  Description 

■  -  ■  .n  —   ■■!■. 

ALPHA  Assumed  junction  angle  of  0°,  degree 

DIAM  Duct  diameter,  foot 

FRCTL  Friction  loss  coefficient  for  athwartships 

length 

KOUNT  Working  variable  for  craft  condition  I 

AJCT  Length  of  fore  and  aft  ducting  to  pump  inlet, 

foot 

FRCTJ  Friction  loss  coefficient  for  junction 

XLAMD  Working  variable  for  AMIXL  calculation 

FLONG  Length  of  athwartships  ducting,  foot 

AJCTL  Total  less  coefficient  of  junction 

AMIXL  Mixing  loss  coefficient  of  junction 

V  Average  velocity,  foot  per  second 

XPUMP  Number  of  gas  turbines  +  1 

RE  Reynolds  number  of  duct 

BETA (J)  Data  array  of  junction  angles,  degree 

ALAMD(J)       Data  array  of  mixing  loss  coefficient 

corresponding  to  BETA (J) 

FRICT  Function  statement,  calculates  Moody  friction 

factor  for  smooth  pipe,  based  on  Reynolds 
number 
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LIST  OP  VARIABLES  USED  IN  PIPE 
NOT  IN  COMMON 

Variable  Description 

KOUNT  Working  variable  for  craft  condition  I 

V  Average  velocity,  foot  per  second 

DIAM  Duct  diameter,  foot 

RE  Reynolds  number  of  duct 

APIPE  Length  of  fore  and  aft  pipe  to  pump  inlet, 

foot 

FRCTL  Friction  loss  coefficient  for  athwartships 

length 

XPUMP  Number  of  gas  turbines 

XKT  Friction  loss  coefficient  of  fore  and  aft 

length 

XLONG  Length  of  athwartships  pipe,  foot 

FRICT  Function  statement,  calculates  Moody  friction 

factor  for  smooth  pipe,  based  on  Reynolds 
number 
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LIST  OF  VARIABLES  USED  IN  DIVRG 
NOT  IN  COMMON 

Variable  Descrintion 

.1        i  i    i    i  *■     i  i 

FRICT  Function  statement,  calculates  Moody  friction 

factor  for  smooth  pipe,  based  on  Reynolds 
number 

KOUNT  Working  variable  for  craft  condition  I 

RE  Reynolds  number  of  duct 

DIVL  Divergence  loss  coefficient,  not  including 

friction 

DWGT1  Duct  weight  of  divergence  angle,  pound 

DWGT2  Duct  weight  of  divergence  length,  pound 

DWGT3  Duct  weight  of  pump  inlet  angle,  pound 

WWGTl  Water  weight  of  divergence  angle,  pound 

WWGT2  Water  weight  of  divergence  length,  pound 

WWGT3  Water  weight  of  pump  inlet  angle,  pound 

CGWX1  Longitudinal  center  of  gravity  of  divergence 

angle,  from  transom,  foot 

CGWX2  Longitudinal  center  of  gravity  of  divergence 

length,  from  transom,  foot 

CGWX3  Longitudinal  center  of  gravity  of  pump  inlet 

angle,  from  transom,  foot 

DIAM  Duct  diameter,  foot 

FRCTL  Friction  loss  coefficient  for  athwartships 

length 

DIVLC  Total  divergence  loss  coefficient 

XPUMP  Number  of  gas  turbines  +  1 

ANGLE  Divergence  angle,  also  pump  inlet  angle,  radian 

HEADL  Total  divergence  head  loss,  foot 

FLONG  Athwartships  length,  foot 

-99- 


Variable 

XLONG 

V 

ADIV 

THETA(J) 

COEF(J) 


Description 

Fore  and  aft  length,  foot 

Average  velocity,  foot  per  second 

Divergence  length,  ADIV  =  X LONG/COS (ANGLE) , foot 

Data  array  of  divergence  angles,  degree 

Data  array  of  divergence  loss  coefficients 
without  friction,  corresponding  to  THETA(J) 
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LIST  OF  VARIABLES  USED  IN  NOZZL 
NOT  IN  COMMON 

Variable  Description 

KOUNT  Working  variable  for  craft  condition  I 

XLPS  Dummy  nozzle  length,  foot 

XFAC  Check  on  whether  nozzle  exits  through  bottom 

or  stern 

XPUMP  Number  of  pumps 

ANOZ  Optimum  nozzle  depression  angle,  radian 

XLNOZ  Nozzle  length,  foot 

DT  Nozzle  throat  diameter,  foot 

DJ  Nozzle  jet  diameter,  foot 

XCORR  Nozzle  head  loss  factor 

RE  Reynolds  number,  based  on  average  diameter 

and  velocity 

AREAl  Throat  pipe  area,  square  foot 

AJET1  Jet  pipe  area,  square  foot 

QQ  Flow  rate  per  nozzle,  cubic  foot  per  second 

FRICT  Function  statement,  calculates  Moody  friction 

factor  for  smooth  pipe,  based  on  Reynolds 
number 
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LIST  OF  VARIABLES  USED  IN  PUMP 
NOT  IN  COMMON 

Variable  DescriDtion 

1  —     ■    M  I         I  ■-■■■!■■  -      -  -  *~ 

K  Cruise   condition  indicator 

J  Take-off   condition   indicator 

I  Index    for  do   loops 

XPUMP  Number  of  pumps 

THOMI  Lower   limit   of  Thoma's   criterion    for   single 

inducer   axial  pump 

DRAT  Hub   to  tip   diameter  ratio 

XXLP  Factor   for  pump    length 

CW  Weight   coefficient 

QX  Ratio   of   cruise    flow   rate   to  take-off    flow 

rate 

CA  Factor   for   RPM  calculation 

CB  Factor   for  RPM  calculation 

HX  Patio  of  design  pump   head  to  off  design  head 

RX  Ratio  of   design   to   off   design  pump    RPM 

PHRAT  Ratio  of  off  design   to   design    flow 

coefficients 

XNGT  Number  of  gas    turbines 

ETAPP  Product   of  pump   and  gearbox   efficiencies 

HP  Inducer  head   for  one    and   two   stage    axial 

pump   designs,    foot 

THOMS  Thoma's    cavitation    criterion    for   inlet   to 

axial   stage 

HHP  Axial  stage  pump  head,  foot 

BETA2  Exit  blade  angle,  radian 
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Variable 


Description 


XNNS 

BD 

CC 

AA 

IMPL 

M 

N 

PC(K,J) 

PCA(K,J) 

PCC(K,J) 

XRPM(M) 

XD1(M) 

XD2 (M) 

RPK(M) 

XPUP(M) 

YLP  (M) 

XERAT(M) 

APUP(NSTG) 
WRAT(NSTG) 


Non-dimensional   specific   speed 

Impeller  exit  width   ratio 

Factor  in    flow   coefficient   calculation 

Factor  in    flow   coefficient   calculation 

Maximum  number  of   impellers  permitted   for 
centrifugal  pump,    =   10 

Working  variable    for  number  of   centrifugal 
pump   impellers 

Working  variable    for  NSTG 

Inducer  head  curve    coefficients 

Inducer  plus   axial   stage  head  curve    coefficients 

Centrifugal   pump  head  curve   coefficients 

Working  variable    for  off  design   PPM   of 
centrifugal  pump  with   M  impellers,    RPM 

Working  variable    for   inlet  tip   diameter,    D1S, 
foot 

Working  variable    for  exit  tip   diameter,    D2S, 
foot 

Working   variable    for   design   RPM  of   centrifugal 
pump  with  M  impellers,    square    foot 

Working  variable    for   centrifugal  pump   area 
with  M  impellers,    foot  squared 

Working  variable    for   centrifugal  pump   length 
with  M  impellers,    foot 

Working  variable    for  gear   ratio   required   for 
centrifugal  pump  with  M  impellers 

Inlet   area  of  pump   NSTG,    square    foot 

Weight   ratio  of  pump   NSTG,    including  pump   dry 
and  wet  weight,    gearbox   and   fuel 
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Variable  Description 

WD(M)  Working  variable  for  XWD(NSTG)  for 

centrifugal  pump  with  M  impellers 

WW(M)  Working  variable  for  XWW(NSTG)  for  centri- 

fugal pump  with  M  impellers 

WWG<M)         Working  variable  for  XWG(NSTG)  for  centri- 
fugal pump  with  M  impellers 
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LIST  OF  VARIABLES  USED  IN  FUEL 
NOT  IN  COMMON 

Variables  Description 

CFS  Constant  for  SFC  calculation 

CA  Cruise  condition  drag  to  lift  ratio 

CD  1  +  total  system  head  loss  coefficient, 

based  on  jet  velocity 

N  Number  of  intervals  endurance  is  divided  into 

+  1 

XN  Number  of  intervals  endurance  is  divided  into 

TI  Time  to  cover  one  ranae  interval  at  constant 

V0(1),  hour 

IJK  Index  for  shifting  SFC  curves 

XJ  Factor  to  convert  SFC  if  SHP  is  less  than  70% 

of  design  SHP 

I  Working  index  for  N 

H  Head  required,  foot 

SHPP  Total  thrust  required,  horsepower 

SHNG  Total  thrust  required  per  engine,  horsepower 

M  Working  index  for  N 

WT(ITIME)      Weight  of  fuel  used  in  time  increment  ITIME, 

pound 

DIS(ITIME)  Displacement   at  time    increment   ITIME,    pound 

VJJ(ITIME)  Jet  velocity   at  time   increment   ITIME,    foot 

per  second 
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LIST  OF  VARIABLES  USED  IN  PTTRN 
NOT  IN  COMMON 

Variable  Description 

PSI  (N)         Current  basepoint  coordinate  of  parameter  N 

THETA(N)       Previous  basepoint  coordinate  of  parameter  N 

PHI (N)         Present  exploratory  point  coordinate  of 

parameter  N 

DEL(N)  Current  step  size  of  parameter  N 

DELMIN(N)  Minimum  step  size  of  parameter  N 

DIR(N)  Last  successful  direction  of  parameter  N 

SAVE(N)  Working  variable  for  PHI(N) 

S  Working  variable  for  function  value 

SPHI  Working  variable  for  function  value  at  PHI 

coordinates 

SPSI  Current  best  function  value  at  PSI  coor- 

dinates 

NUMB  Counter  for  minimum  step  size  check 

RHO  Step  size  change  factor 

ICALL  Indicator  of  current  point  move 

K  Index  for  do  loop 

I  Index  for  do  loop 

N  Number  of  parameters  of  search 

SIGN(N)  Directed  step  size  of  parameter  N 
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LIST  OF  VARIABLES  USED  IN  OUTPUT 
NOT  IN  COMMON 


Variable 

IPRNT 

J 

SUM 

I 

DRATO 

IK 

XNGT 

K 

L 

LL 

KLL 

M 

NIMP 

TFM 

WTRAT 

WRATF 

ENGN(IENGN) 

VJRAT(I) 

VI RAT (I) 


Description 

Print  data  set  reference  number 

Working  variable  for  condition  I 

Total  duct  head  loss,  excluding  elevation, 
foot 

Index  for  do  loop 

Strut  diffuser  area  ratio 

Index  for  number  of  elbows  in  system 

Number  of  gas  turbines 

Index  for  do  loop 

Index  for  implied  do  loop  in  output  statement 

Working  variable  to  point  to  correct  head 
loss  for  output 

Working  variable  to  point  to  correct  format 
statement  for  head  loss  output 

Index  for  implied  do  loop  in  output  statement 

Number  of  impellers  in  pump 

Strut  thickness  at  flying  waterline,  foot 

Total  propulsion  system  weight  ratio 

Propulsion  system  weight  ratio,  excluding  fuel 

Contains  name  of  engine  IENGN 

Jet  velocity  ratio  at  craft  condition  I 

Inlet  velocitv  ratio  at  craft  condition  I 


HEADL(I,ICOMP)  Head  loss  of  component  ICOMP  at  craft  con- 
dition I 


C0NDS(2,I) 


Label  for  craft  condition  I 
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Variable  Description 

ELBWS(2,IK)  Label  for  elbow  IK 

PC (I)  Propulsive  coefficient  at  craft  condition  I 

LABEL (5 ,M)  Labels  for  output 
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